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The role of the magnetoelectric effect upon optical reflectivity is studied by adapting an 
electrodynamic-based model for a system composed by a 2D metallic film in contact with an 
extended multiferroic material exhibiting weak ferromagnet ism. The well-known Nakayama's 
boundary condition is reformulated by taking into account the magnetoelectric coupling as well 
as an externally applied magnetic field B in arbitrary direction. It is found that the relative 
reflectance shows strong fluctuations for in-plane B -fields and might be controllable under its 
azimuthal rotation in the THz regime, in particular, close to the antiferromagnetic resonance 
frequency associated to the multiferroic slab. 
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I. INTRODUCTION 

Magnetoelectric (ME) effects in multiferroic (MF) or 
ferromagnetic (metallic) films have brought remarkable 
interest since promising technological applications in 
spintronics and ultrafast electric field control on mag- 
netic data storage are seen as imminentPp. Charac- 
terization of the relative strength for the ME coupling 
can be obtained by implementing terahertz spectroscopy 
in rare earth manganites of the type RMn03 (R=Tb, 
Gd, Dy, Eu:Y]PPPp demonstrating that the generated 
electromagnons (mixed spin-waves and photon states) 
represent, among others, the signature of the ME ef- 
fect for an approximate range of frequencies between 
10 cm -1 to 40 cm -1 at temperatures where antiferro- 
magnetic resonance modes (AFMR) coexist, or more 
recently, the key mechanism for controllable magne- 
tochromism in Ba2Mg2Fei2 022 hexaferrites^. The pri- 
mary origin for the ME coupling is commonly associ- 
ated to the Dzyaloshinskii-Moriya relativistic exchange- 
interact iorPP which is appropriate for the description 
for asymmetric spin wave dispersion on double layer 
Fe-films^ as well as for those materials where weak 
ferromagnetism emerges, namely the Ilmenite FeTiOs, 
TbMn0 3 , Eui-^MnOs (0 < x < 0.3 at T < 40 KP) or 
the widely studied pyroelectric ferromagnet BaMnFj^. 
Weak ferromagnetism on this compound is generated 
by canting effects between antiferromagnetic sub-lattices, 
leading into a spontaneous polarization P perpendicular 
to the resulting magnetization M 13 . Considerations in 
the symmetry change of the static polarization and mag- 
netization fields have brought interesting unconventional 
optical phenomena labeled as non-reciprocal dichroism 
associated to the sign reversal of P x M, recently re- 
ported in the perovskite Euo.55Yo.45Mn03, with magne- 
toelectric activity for photon energies around 0.8 meV 



(sub THz regime) in the cycloidal phase at 4 KP=l In- 
tense activity in the last decade has also been dedicated 
to achieve possible optical and photonic band gap con- 
trol via Surface Plasmon (SP) propagation in periodic 
arrays 15 , since modern lithographic techniques allows to 
design functional objects with almost any geometrical 
pattern at sub- wavelength scale 16 . Plasmon localization 
and its coupling with incident light depend on the di- 
electric properties of the metal in conjunction of its sur- 
rounding environment, enlightening an alternative route 
for engineering highly efficient SP photonic devices via 
externally applied fields, rare earth doping or electron 
charge transference from the modified metaP^. In this 
communication we study an electrodynamic-based model 
for estimating the optical response generated by the con- 
tact between a material exhibiting weak ferromagnetism 
in contact with a 2D metallic film. It is found that a spe- 
cific strength of the ME interaction might couple with 
localized charge-sheet modes for electron carrier densi- 
ties close to 10 16 cm -2 and incident frequencies around 
18 cm -1 , leading into a change in the reflectance from the 
metallic film. Applied magnetic field effects on relative 
reflective are also discussed in section III. 

II. MODEL 

Localized charge-sheet modes in a 2D conducting 
medium in the framework of Drude appro ximation is ob- 
tained from the Nakayama resurp 18 f 19 f 2Q t 
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where Kj = (qy — SjUO 2 /c 2 ) , (j = 1,2) corresponds 
to the quasiwavevector in 7- direction, Q$ is defined 
as ve 2 /some and v denotes the electron density con- 
centration in a two dimensional space. The term Sj 
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represents the relative dielectric function value for j-th 
medium, with e\ = 1 for vacuum. In the the range 
of wavelengths behind the far infrared radiation (< 1 
mm), the dielectric function approaches to the well rec- 
ognized Lyddane- Sachs-Teller (LST) relationship: e 2 ~ 
(1 + Xoo) (<^lA^t) 2 , where Xoo corresponds to the dielec- 
tric permittivity of the medium j = 2 and Wl,(t) re P _ 
resents the longitudinal (transverse)-optical phonon fre- 
quency. For numerical purposes, we have set (uol/ojt) 2 ~ 
1.07, which coincides with the relationship for the 6-axis 
normal phonon modes in BaMnF4. The permittivity 
Xoo is a functional depending on mechanical strain de- 
formations and polarization field depletion in the prox- 
imities between the multiferroic slab and metal film^J 
and is taken as constant for zero applied (electric) field 
and fixed temperature. Formula ([!]) is derived by solv- 
ing the complete set of Maxwell equations with normal 
(TM wave) incidence for Z > 0, and boundary condi- 
tions on the plane Z = with the ansatz for prop- 
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in the region 



Z < 0. Magnetoelectric effects are taken into consid- 
eration throughout the transverse susceptibility x me and 
the electric displacement vector D is written into the 
constitutive equation like D = + 47rx me H. Af- 
ter inserting the additional term 47rx me H, the expres- 
sion shall be modified under k 2 —> &2 + 47nc<jx me / c - 
In the plane Z = 0, and in agreement with the geo- 
metrical configuration showed in Figure ([!]), the non- 
zero surface current density component is defined as 
Jy = crEy, where a corresponds to the ayy-element 
of the generalized conductivity tensor^, and Ey is the 
electrical field propagating on Z direction. The generic 
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FIG. 1: Conducting Charge-Sheet in contact with a multifer- 
roic surface. The polarization vector P and the wavevector 
of coupled excitations qy are also depicted in the diagram. 
Weak ferromagnetic magnetization vector M is produced by 
interacting antiferromagnetic sublattices with relative canting 
angle 6c- 

expression for the transverse susceptibility x me is ob- 
tained from first principles 23 , 24 and it can be summarized 

as: 4:7riux r ' 

where g = 
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eter which is an involved function of the canting an- 
gle between two adjacent (antiferromagnetic) sublattices, 
the spontaneous magnetization M and the polarization 
vector P, as well as the parameters (jJ m ( p ). g is con- 
veniently given in units of mm -2 all throughout this 
paper 25 , in concordance with the spectral weight intrin- 
sically associated to the fitting procedure for the trans- 
mittance spectra via Lorentzian model in various multi- 
ferroic species, namely RMn 2 5 (R:Y,Tb), TbMn0 3 or 
LuMnOs 26 , and its dependence with the externally ap- 
plied magnetic field has been neglected for small canting 
angles (See for instance Eqns. (38) and (47) in ref. (P*J)). 
Two main poles are clearly identified for x me : the opti- 
cal antiferromagnetic resonance mode (AFMR) uj m and 
the soft-phonon along M with resonance frequency uo pi 
with uj p > uo m . Classical plasmon excitations in low 2D 
carrier electron density are experimentally detected and 
theoretically estimated for wavevectors q < 1.4 cm -1 and 
frequencies uj < 0.5 meV 27 f 28 f 29 l, therefore the condition 
q Y » SjUJ 2 /c 2 remains valid in the range of interest, 
and the dispersion relationship for the coupled magneto- 
electric plasma mode is obtained by solving the modified 
equation ([!]): 
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with Q = 47rzco>x me / c + 71 (uj/2ttc) , 71 = 
47r 2 c(e 1 +e 2 )/n s and 72 = 167r 2 ce 1 /n s . For X me = 0, 
i.e., no magnetoelectric effects taken under consider- 
ation, we reproduce the expression for the localized 
plasmon moc 
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g (0c, M, P, co> m , uj p ) is a coupling param- 



where (+) sign in equation has been selected. Com- 
plex index of refraction ft (a;) is directly estimated from 
the wavenumber 3 -^ qy : h(uj) = cqy/uj. The lowest-order 
reflectance coefficient R (uj) for normal incidence is de- 
fined as R(uj) =| n(uj) — 1 | 2 / | n(uj) + 1 | 2 and its 
numerical profile discussed on the next section. 

III. RESULTS AND DISCUSSION 

Figure Q exhibits the reflectance response as a func- 
tion of the 2D electronic carrier concentration for 
different wavelengths and the magneto-electric coupling 
parameter g fixed at 0.02. Dotted curve (a) is set as 
reference for g = 0. In this case, reflectivity minima 
are located at wavelengths A c = 2tt (si + £2) c/fis, or 
A" 1 oc indicating that the critical wavelength for bare 
plasmon excitations is larger as the electronic concen- 
tration decreases. The AFMR mode lies in the range 
THz range, with uj m ~ 0.54 THz, while the transverse 
phonon frequency is taken as 7.53 THz for the BaMnF4 
compound 31 . Metallic behavior predominates for con- 
centrations higher than 10 18 cm -2 and smaller than 10 16 
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cm -2 and selected wavelengths between 0.4 mm and 0.7 
mm. Plasmonic resonance effects are important for car- 
rier densities around 10 1T cm -2 , where radiative absorp- 
tion or antireflective phenomena become strong and the 
reflectance spectrum is therefore significantly modified 
by diminishing the percentage of absorbed radiation only 
when the external frequency approaches the characteris- 
tic mode o; m , and g ^ 0. Figure (|3| depicts the shift- 
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FIG. 2: Zero field reflectance response as a function of elec- 
tron carrier density v. Curves (a) and (b) correspond to re- 
sponses with g = and g ~ 0.02 respectively, close to the 
optic antiferromagnetic resonance mode cj m . 

ing of the minimum of reflectance in the (V, g) plane. 
The ME effect becomes relevant by decreasing the criti- 
cal carrier density v c as g increases, and it remains un- 
modified for those frequencies far away from the AFMR 
characteristic mode as indicated in lines (d) and (e). v c 
shall be understood as the electron carrier concentration 
which maximizes the antireflective effect for the compos- 
ite metal/multiferroic system. Applied magnetic field B 
enters into the formalism by taking symmetry considera- 
tions upon the dependence of the electrical conductivity 
as a function of B under the transformation a — >• a (£?), 

with a (B) = iQs c ~ luj (<^ 2 — <^#) > anc ^ °°b = eB/m as 
the cyclotron frequency. Expression Q changes as qp = 
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with Q' = Q — jiLOg/J^Trc) 2 . The classical localized mag- 



netoplasmon mode 
Bafl 



(27TCT 

(pi) is rewritten for g 

4:7T 2 C 2 qY 

7i 



(" 2 -<4) 
(2ttc) 2 



(4) 
ag- 

and under 




(5) 



in similarity with result In this particular case, 

the antireflective condition (n = 1) depends on the 
external magnetic field intensity as A" 1 = 71-/71 + 

\/(7r/7i) 2 + (cjb/27tc) 2 , which leads into a quadratic 



correlation A" 1 ex B 2 for j 1 ujb/2tt 2 c << 1. For 
an orientation of B, equation shall be modi- 
fied on its right side accordingly with Qs CUJ ~ 2 — » 
Qs c (u 2 — ojg) F (rix, ny , where F (•) is a func- 
tion of the directors nx,Y,z 33 - Maximum change in the 
relative reflectivity takes place for B \\ X-axis (Voigt 
configuration) and parallel to M, while no significant 
(or zero) variation in AR/R = R (B) /R(0) — l 34 comes 
out for B || F-axis, suggesting that an effective rotation 
on the weak magnetization vector M generates quasi- 
periodic fluctuations upon the optical activity close to 
the antiferromagnetic resonance frequency and dimin- 
ishes coherence in the collective excitations propagating 
along Y- direction, as shown in figure Q. In-plane ap- 
plied field B effects on the reflectance as a function of 
carrier density v are illustrated in Fig. ([5]). R (B) tends 
to augment for B parallel to X-axis and decreases for B 
along —X axis. Curve (b) for null B overlaps the out- 
come of R for B = 1.5 T, <ft = tt/2 (i.e., parallel to Y 
axis) and the same g parameter, indicating no substan- 
tial variation in the optical reflectance for applied fields 
in the same direction of the plasmonic wavevector qy. 
Strong fluctuations around the baseline AR = in the 
relative reflectance factor AR/R are mainly observed for 
wavelengths close to the AFMR frequency 27rc/cj m , and 
their intensities tend to decrease for stronger ME cou- 
plings (See figures ([6| and (|8|). These effect does not 
appear when the g is set as zero (not shown) , suggesting 
a direct correlation between the dressed confined plas- 
monic excitations at the Drude metal, the intensity of 
the applied field and its optical response. Dotted line (e) 
in Fig. ^ is set as reference for B perpendicular to the 
plane Z = for parameters corresponding to curve (b) in 
Figure (|6|. Figure ^ shows the reflectivity response for 
externally applied fields parallel to the Z — plane as a 
function of the incident wavelengths. Iso-reflective lines 
for AR/R close to 2i\c/uj rn and the externally applied 
magnetic field (in Z direction) are shown in Figure (|9|. 
Projected lines preserve symmetrical distribution under 
magnetic field inversion but it changes for wavelengths 
around 27Tc/uj m and g > 0, i.e., regardless the applied 
magnetic field orientation, the relative reflectance flips its 
sign at 2i\c/uj rn . The reflectance function R (B) reaches 
values from 0.2R(B = 0) to 0.8R(B = 0) for intensities 
lying in the interval 0.52 < B < 2.3T. On the edge of the 
THz bandwidth, i.e. A# ~ 1 mm, the reflectance factor 
R(Xe) remains constant (around 0.1) for magnetic field 
intensities lesser than 3 T and augments exponentially 
up to - 0.4 for 3 < B < 8 T and g = 0.02. In the limit 
of the microwave band (Am °° 10 mm), R(Xm) changes 
from 0.8 to 0.96 for < B < 1 T with no significant 
influence of ME coupling on the metallic behavior. 

IV. CONCLUDING REMARKS 

We have developed a model for studying the magneto- 
electric interactions on 2D plasmonic modes in the THz 
range for a metal/multiferroic device. The multiferroic 
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FIG. 3: Critical carrier density v c as a function of the ME 
coupling parameter g for different wavelengths. v c is strongly 
depending on g only for external frequencies close to AFMR 
mode ujm- 



numerical values were taken for BaMnF4, showing 
that a particular condition for antireflectivity might be 
adjustable by varying the intensity of the applied field, 
its orientability or the incident frequency of radiation in 
a range lj > cj m , while the zero field reflectance increases 
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60% for a change in the ME coupling around 
= 0.02 mm" 2 . The quantity R(B) at B = IT 
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FIG. 4: Relative Reflectance factor for in-plane applied mag- 
netic field close to AFMR frequency at 2tyc/uj = 0.53 mm, at 
g = 0.02 as a function of azimuthal angle (f). 



FIG. 5: Reflectance as a function of electron carrier density 
for in-plane applied field B and 2ttc/u = 0.53 mm. Line (a) 
corresponds to uncoupled (g = 0) response. 



parallel to the X axis reaches values up to 10R (0) in 
contrast with approximately 1.2^(0) for the same field 
intensity but for B perpendicularly applied to the XY 
plane, indicating notorious enhancement in the value 
of the relative reflectance in the former case. Although 
all numerical simulations were conducted for S2 — Ke\ 
(K being a positive number), simultaneous electric field 
control Eo on optical properties for the composite device 
might also be achieved under the dielectric function 
dependence for a multiferroic material £2p(Eo)], the 
polarization P (Eo) and temperature, issue that shall be 
addressed in further investigations. 



medium exhibits weak ferromagnetism and the metallic 
behavior enters into the formalism in the framework of 
the classical Drude-Lorentz model. Relative reflectance 
response for normal incidence is numerically calculated 
for different ME coupling strengths and wavelengths 
close to the optical antiferromagnetic resonance fre- 
quency uj m . Characteristic soft phonon and AFMR 
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